Mathematical modeling of the metastatic process

Mathematical modeling in cancer has been growing in popularity and impact since its inception in 1932. The first theoretical mathematical modeling in cancer research was focused on understanding tumor growth laws and has grown to include the competition between healthy and normal tissue, carcinogenesis, therapy and metastasis. It is the latter topic, metastasis, on which we will focus this short review, specifically discussing various computational and mathematical models of different portions of the metastatic process, including: the emergence of the metastatic phenotype, the timing and size distribution of metastases, the factors that influence the dormancy of micrometastases and patterns of spread from a given primary tumor.Comment: 24 pages, 6 figures, Revie

Modelling of Tirapazamine effects on solid tumour morphology

Bioreductive drugs are in clinical practice to exploit the resistance from tumour microenvironments especially in the hypoxic region of tumour. We pre-sented a tumour treatment model to capture the pharmacology of one of the most prominent bioreductive drugs, Tirapazamine (TPZ) which is in clinical trials I and II. We calculated solid tumour mass in our previous work and then integrated that model with TPZ infusion. We calculated TPZ cytotoxicity, concentration, penetra-tion with increasing distance from blood vessel and offered resistance from micro-environments for drug penetration inside the tumour while considering each cell as an individual entity. The impact of these factors on tumour morphology is also showed to see the drug behaviour inside animals/humans tumours. We maintained the heterogeneity factors in presented model as observed in real tumour mass es-pecially in terms of cells proliferation, cell movement, extracellular matrix (ECM) interaction, and the gradients of partial oxygen pressure (pO2) inside tumour cells during the whole growth and treatment activity. The results suggest that TPZ high concentration in combination with chemotherapy should be given to get maximum abnormal cell killing. This model can be a good choice for oncologists and re-searchers to explore more about TPZ action inside solid tumour

Using DNA Methylation Patterns to Infer Tumor Ancestry

Background: Exactly how human tumors grow is uncertain because serial observations are impractical. One approach to reconstruct the histories of individual human cancers is to analyze the current genomic variation between its cells. The greater the variations, on average, the greater the time since the last clonal evolution cycle (‘‘a molecular clock hypothesis’’). Here we analyze passenger DNA methylation patterns from opposite sides of 12 primary human colorectal cancers (CRCs) to evaluate whether the variation (pairwise distances between epialleles) is consistent with a single clonal expansion after transformation. Methodology/Principal Findings: Data from 12 primary CRCs are compared to epigenomic data simulated under a single clonal expansion for a variety of possible growth scenarios. We find that for many different growth rates, a single clonal expansion can explain the population variation in 11 out of 12 CRCs. In eight CRCs, the cells from different glands are all equally distantly related, and cells sampled from the same tumor half appear no more closely related than cells sampled from opposite tumor halves. In these tumors, growth appears consistent with a single ‘‘symmetric’ ’ clonal expansion. In three CRCs, the variation in epigenetic distances was different between sides, but this asymmetry could be explained by a single clonal expansion with one region of a tumor having undergone more cell division than the other. The variation in one CRC was complex and inconsistent with a simple single clonal expansion

Biological and clinical significance of cancer stem cell plasticity

In the past decade, the traditional view of cancers as a homogeneous collection of malignant cells is being replaced by a model of ever increasing complexity suggesting that cancers are complex tissues composed of multiple cell types. This complex model of tumorigenesis has been well supported by a growing body of evidence indicating that most cancers including those derived from blood and solid tissues display a hierarchical organization of tumor cells with phenotypic and functional heterogeneity and at the apex of this hierarchy are cells capable of self-renewal. These “tumor imitating cells” or “cancer stem cells” drive tumorigenesis and contribute to metastasis, treatment resistance and tumor relapse. Although tumor stem cells themselves may display both genetic and phenotypic heterogeneity, recent studies have demonstrated that cancer stem cells maintain plasticity to transition between mesenchymal-like (EMT) and epithelial-like (MET) states, which may be regulated by the tumor microenvironment. These stem cell state transitions may play a fundamental role in tumor progression and treatment resistance. In this review, we discuss the emerging knowledge regarding the plasticity of cancer stem cells with an emphasis on the signaling pathways and noncoding RNAs including microRNAs (miRNA) and long non-coding RNAs (lncRNAs) in regulation of this plasticity during tumor growth and metastasis. Lastly, we point out the importance of targeting both the EMT and MET states of CSCs in order to eliminate these lethal seeds of cancers. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1186/s40169-014-0032-3) contains supplementary material, which is available to authorized users

The role of Allee effect in modelling post resection recurrence of glioblastoma

Resection of the bulk of a tumour often cannot eliminate all cancer cells, due to their infiltration into the surrounding healthy tissue. This may lead to recurrence of the tumour at a later time. We use a reaction-diffusion equation based model of tumour growth to investigate how the invasion front is delayed by resection, and how this depends on the density and behaviour of the remaining cancer cells. We show that the delay time is highly sensitive to qualitative details of the proliferation dynamics of the cancer cell population. The typically assumed logistic type proliferation leads to unrealistic results, predicting immediate recurrence. We find that in glioblastoma cell cultures the cell proliferation rate is an increasing function of the density at small cell densities. Our analysis suggests that cooperative behaviour of cancer cells, analogous to the Allee effect in ecology, can play a critical role in determining the time until tumour recurrence

Bcl-2 inhibits apoptosis by increasing the time-to-death and intrinsic cell-to-cell variations in the mitochondrial pathway of cell death

BH3 mimetics have been proposed as new anticancer therapeutics. They target anti-apoptotic Bcl-2 proteins, up-regulation of which has been implicated in the resistance of many cancer cells, particularly leukemia and lymphoma cells, to apoptosis. Using probabilistic computational modeling of the mitochondrial pathway of apoptosis, verified by single-cell experimental observations, we develop a model of Bcl-2 inhibition of apoptosis. Our results clarify how Bcl-2 imparts its anti-apoptotic role by increasing the time-to-death and cell-to-cell variability. We also show that although the commitment to death is highly impacted by differences in protein levels at the time of stimulation, inherent stochastic fluctuations in apoptotic signaling are sufficient to induce cell-to-cell variability and to allow single cells to escape death. This study suggests that intrinsic cell-to-cell stochastic variability in apoptotic signaling is sufficient to cause fractional killing of cancer cells after exposure to BH3 mimetics. This is an unanticipated facet of cancer chemoresistance.Comment: 11 pages, In pres

Number and Size Distribution of Colorectal Adenomas under the Multistage Clonal Expansion Model of Cancer

Colorectal cancer (CRC) is believed to arise from mutant stem cells in colonic crypts that undergo a well-characterized progression involving benign adenoma, the precursor to invasive carcinoma. Although a number of (epi)genetic events have been identified as drivers of this process, little is known about the dynamics involved in the stage-wise progression from the first appearance of an adenoma to its ultimate conversion to malignant cancer. By the time adenomas become endoscopically detectable (i.e., are in the range of 1–2 mm in diameter), adenomas are already comprised of hundreds of thousands of cells and may have been in existence for several years if not decades. Thus, a large fraction of adenomas may actually remain undetected during endoscopic screening and, at least in principle, could give rise to cancer before they are detected. It is therefore of importance to establish what fraction of adenomas is detectable, both as a function of when the colon is screened for neoplasia and as a function of the achievable detection limit. To this end, we have derived mathematical expressions for the detectable adenoma number and size distributions based on a recently developed stochastic model of CRC. Our results and illustrations using these expressions suggest (1) that screening efficacy is critically dependent on the detection threshold and implicit knowledge of the relevant stem cell fraction in adenomas, (2) that a large fraction of non-extinct adenomas remains likely undetected assuming plausible detection thresholds and cell division rates, and (3), under a realistic description of adenoma initiation, growth and progression to CRC, the empirical prevalence of adenomas is likely inflated with lesions that are not on the pathway to cancer

Quantitative Interpretation of a Genetic Model of Carcinogenesis Using Computer Simulations

The genetic model of tumorigenesis by Vogelstein et al. (V theory) and the molecular definition of cancer hallmarks by Hanahan and Weinberg (W theory) represent two of the most comprehensive and systemic understandings of cancer. Here, we develop a mathematical model that quantitatively interprets these seminal cancer theories, starting from a set of equations describing the short life cycle of an individual cell in uterine epithelium during tissue regeneration. The process of malignant transformation of an individual cell is followed and the tissue (or tumor) is described as a composite of individual cells in order to quantitatively account for intra-tumor heterogeneity. Our model describes normal tissue regeneration, malignant transformation, cancer incidence including dormant/transient tumors, and tumor evolution. Further, a novel mechanism for the initiation of metastasis resulting from substantial cell death is proposed. Finally, model simulations suggest two different mechanisms of metastatic inefficiency for aggressive and less aggressive cancer cells. Our work suggests that cellular de-differentiation is one major oncogenic pathway, a hypothesis based on a numerical description of a cell's differentiation status that can effectively and mathematically interpret some major concepts in V/W theories such as progressive transformation of normal cells, tumor evolution, and cancer hallmarks. Our model is a mathematical interpretation of cancer phenotypes that complements the well developed V/W theories based upon description of causal biological and molecular events. It is possible that further developments incorporating patient- and tissue-specific variables may build an even more comprehensive model to explain clinical observations and provide some novel insights for understanding cancer